Optical conductivity due to orbital polarons in systems with orbital degeneracy 
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(i) are determined by the string picture, i.e., flipped orbitals along the hole hopping path, and 
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the optical transitions. 
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I. INTRODUCTION 

Electronic excitations at finite frequency are a com- 
mon feature of doped Mott or charge-transfer insula- 
tors. Such excitations were observed experimentally in 
high temperature superconductors by optical absorption 
<j(lo) shortly after these systems were discovered^ and 
were extensively studied in the t-J model by several 
groups, 2-6 providing valuable insights into the charge dy- 
namics in doped cuprates. It was also realized that con- 
siderable transfers of spectral weight in the optical spec- 
tra occur and new states arise which generate spectral 
intensity within the Mott-Hubbard gap^— These stud- 
ies have shown that the key assumption of the Fermi 
liquid theory that the low-energy excitation spectrum 
stands in a one-to-one correspondence with that of a 
non-interacting system has to be revised when the elec- 
trons interact strongly. For instance, drastic deviations 
from the Fermi liquid picture are obtained, in the nor- 
mal state of the copper-oxide high-temperature super- 
conductors, highlighted by a pseudogap, broad spectral 
features, and the resistivity which increases linearly with 
temperature^ It was recently established that the op- 
tical conductivity for Cu02 planes of high temperature 
superconductors exhibits a mid-infrared peak at low dop- 
ing that gradually develops to a band under increasing 
doping and causes an insulator-to-metal transition^ 

Optical conductivity studies play also a very impor- 
tant role in other correlated materials, including Mott 
insulators with active orbital degrees of freedom. These 
systems exhibit rather complex behavior due to the in- 
terrelation between spin, orbital and charge degrees of 
freedom. Recently it was pointed out that the excitations 
to 3z 2 — r 2 orbitals in high-T c cuprates are responsible 
for the observed optical conductivity in the insulating 
stated A study of the undoped three-orbital Hubbard 



model explains the anisotropy of the optical conductiv- 
ity of a pnictide superconductor^ Finally, a well known 
example in this family of compounds are colossal mag- 
netoresitance manganites where their complexity mani- 
fests itself in the large number of competing magnetic 
phases in the phase diagrams of perovskite or layered 
materials^ including the charge-ordered phase at 50% 
doping with the optical conductivity determined by the 
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-like occupied orbital which coexists with charge 



order. 15 In undoped LaMnOs the optical conductivity 
<t(w) shows several features at higher energy, 16 and their 
spectral weights change with temperature. These fea- 
tures and the thermal evolution of their spectral weights 
may be well understood by employing a general relation 
between the spin-orbital superexchange and the spec- 
tral weight distribution in the optical spectra of Mott 
insulators^ This theory is also successful in analyzing 
the temperature dependence of the low-energy optical 
spectra for high spin excitations in LaVC>3,— where spin- 
orbital entanglement plays a role^ 

Detailed investigations of the optical conductivity of 
the ferromagnetic (FM) metallic Lai-aSr^MnOa have 
shown^S (i) a pseudogap in cr(ui) for temperatures above 
the Curie temperature Tq, (ii) the growth of the broad 
incoherent spectrum at low energy < oj < 1.0 cV un- 
der decreasing temperature below Tc, and (iii) a narrow 
Drude peak. In these compounds electronic correlations 
among e g electrons are strong and a orbital liquid stabi- 
lizes the FM metallic phased These features have been 
successfully explained in the theory which focuses on the 
orbital dynamics in a situation when spins do not con- 
tribute and may be neglected^ In this situation orbital 
polarons 2 ^ determine the transport properties and the 
optical spectra. Recently a two-peak structure of the op- 
tical conductivity was discussed for the metallic phase 
of FM manganitesj^i with a far-infrared Drude peak ac- 
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companied by a broad mid-infrared polaron peak. It is 
intriguing whether similar phenomena occur in other or- 
bital systems as well. 

In this paper we investigate the optical conductivity in 
the two-dimensional (2D) orbital model with two active 
tig orbital flavors. This model reveals crucial properties 
we are interested in, and is applicable either to transition 
metal oxides with FM planes and active ti g orbitals when 
the tetragonal crystal field splits off the xy orbital from 
the {yz,zx} doublet filled by one electron at each site, 
as for instance in Sr 2 V04^ or to certain e g planar ma- 
terials such as K2CUF4 or Cs2AgF4 ) 26 ' 27 or to cold-atom 
systems^ with active p orbitals^ Orbital superexchange 
which arises in the strongly correlated regime is responsi- 
ble for alternating orbital (AO) order at half filling. This 
is then the reference state playing a role of the physical 
vacuum below when we consider the optical conductivity 
for a doped Mott insulator. 

The paper is organized as follows. In Sec. [TT| we in- 
troduce the microscopic model and specify typical pa- 
rameters. The model is solved first for a single hole in 
Sec. IIII1 where we analyze the processes of possible hole 
propagation and the role of string states in the optical 
conductivity. Next we present a numerical solution for 
the optical conductivity in Sec. IIV Al and interpret the 
results in Sec. IIVBI A short summary and final conclu- 
sions are presented in Sec. [V] More technical details on 
the performed calculations and on the origin of the op- 
tical transitions that could be of interest only for some 
readers are given in the Appendix „ 

II. GENERALIZED ORBITAL t-J MODEL 

For the purpose of discussing the optical conductiv- 
ity in orbitally degenerate systems we concentrate on 
the recently introduced strong-coupling version of the 
two-orbital Hubbard model for spinlcss fermions on the 
square lattice (when the spins form a FM order and can 
be neglected)^ To be specific, the physical problem to 
which our analysis applies is the FM plane of a Mott insu- 
lator with AO order of ti g (or e g or p) orbitals, i.e., inter- 
acting spinless fermions which undergo one-dimensional 
(ID) nearest neighbor (NN) hopping with conserved or- 
bital flavor: 

H = -t C \,a C 3,a- t C i,b C 3,b + U J2 nianib ■ ^ 

W}||f> W} II a i 

Here cj a and cj b are creation operators for electrons with 
two orbital flavors, and we consider ti g orbitals^ 

\a) = \yz), \b) = \xz), (2) 

labeled by the index of a cubic axis which prohibits the 
electron hopping by symmetry; this notation was in- 
troduced for titanium and vanadium perovskites. 31 The 
summations in Eq. ([T]) are carried over pairs {ij} of NN 
sites in the ab plane. 



In the following we consider an effective Hamiltonian 
obtained by a unitary transformation U which eliminates 
the part of the Hamiltonian that creates/annihilates dou- 
ble occupancies. The purpose is not to eliminate ener- 
getically costly double occupancies, but rather to assess 
accurately to what extent they contribute to the ground 
state energy and wave functions^ The effective Hamil- 
tonian reads 

H t2g = VU^HUV, (3) 

where V is a projection operator that projects the trans- 
formed Hamiltonian on the low energy Hilbert space and 
removes all states with doubly occupied sites. The stan- 
dard perturbation theor y 32 i 33 gives the following expres- 
sion for the generator in Eq. ([3]): 

+ E ( c i,b n ia c j,b ~ c \,b n ™ c i,b)\ + (jji) ■ ( 4 ) 

m\\* V J 

Following this scheme we obtain the generalized orbital 
t-J model with orbital superexchange J = At 2 /U and 
three-site effective next nearest neighbor (NNN) hopping 
t = t 2 /U (both expressions apply when U ^> t)S^ 

H t2g = P{n t + Hj + H% + H ( *>) V , (5) 

Hi = -tJ2 ( C i,6 c i+ft,6 + c i,a c i+ b, Q + H - C ') I ( 6 ) 

1 

n > = \ J Y,( T i T ! -\ n ^i] > ( ? ) 
m v J 

H 3s = ~ T Y1 ( C l-a,fc"i,a c i+a,6 + H - c -) 

i 

- -E( C Ua^ + b, a + H - C -) > ( 8 ) 

i 

i 

~ T E ( C ! T b,a C ^ C lfc C tta,fc + H ' C -) ■ (9) 

i 

Here the summations are carried over sites i € ab plane, 
and the unit vectors {a, b} indicate the bond direction 
in the ab plane. The superexchange J is Ising-like and 
couples NN orbital operators, 

n - \ (n ia - n ib ) , (10) 

on the bonds (ij) in the ab plane. The NN hopping t and 
the effective NNN hopping r contribute only in presence 
of holes as the projection operators V project onto the 
subspace without double occupancies; for more details 
see Ref. HO- As we demonstrate below, a nice feature of 
the Hamiltonian ([SJ is that even on the infinite lattice it 
can be in principle exactly solved in the low energy sector 
by numerical methods. 
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III. A SINGLE HOLE PROBLEM 

We start the analysis of the optical conductivity by 
considering the problem of a single hole. In a Mott 
insulator when there is exactly one electron per lattice 
site only the exchange term Eq. ([7]) contributes and the 
ground state is the 'orbital Neel state' |AO) with AO 
order, playing here a role of the physical vacuum and 
shown schematically in Fig. HJa). Boxes aligned along 
the a (b) lattice direction represent b (a) orbitals, the 
sublattices containing a and b orbitals in the state |AO) 
will be denoted by A and B, respectively. As interactions 
are Ising-like, quantum fluctuations are absent and the 
energy of this state is exactly —\J per bond. This en- 
ergy plays a role of the reference energy of the physical 
vacuum state |AO) in what follows. 

Next we assume that an electron is removed from the 
| AO) state, i.e., a single hole is created in the state with 
AO order, see Fig. [IJb). It will be seen that the motion of 
this hole disrupts the AO order so that the problem has a 
strong similarity with the much studied problem of hole 
motion in an Ising antiferromagnet, sometimes referred 
to as the t-J z model or, more generally, hole motion in an 
antiferromagnet There are two differences: (i) the 
hole motion is directional, i.e., a hole on the A sublatticc 
can move only in b-direction and vice versa, and (ii) the 
term which is usually neglected in the spin t-J model 
but is the only term here responsible for coherent motion 
of the hole. 

We consider the single-hole state c^AO) shown in Fig. 
OJb). In the following we refer to the states shown in Fig. 
[T]by their labels (b), (c), etcetera. Creation of the hole at 
site i raises the expectation value of Hj by AEo = jzJ, 
where z is the number of NNs (here z = 4). The term 
Wt — which in principle has the largest matrix element 
t — couples the state (b) with the state (c). Thereby a 
misaligned orbital is created at i, which further increases 
the expectation value of Hj by AE\ — i(z — 1) J. The 
same holds true for all subsequent hops which involve 
TLt — these create a 'string' of misaligned orbitals, see 
the states (e) and (f). Each of these further misaligned 
orbitals created in step n > 1 increases the energy by 
AE n = |(z — 2) J. The term Tit therefore does not lead 
to a coherent propagation of the hole. Trugman has dis- 
cussed coherent hole motion in the t-J z model whereby a 
hole performs one and a half circular movement along the 
smallest closed loop, i.e., around a plaquette in a square 
lattice.— It is straightforward to see that for this smallest 
2x2 loop this mechanism does not work in the present 
case due to the directional nature of the hopping % t Eq. 
([6]) that excludes the hopping t along closed loops. 

In contrast to this, the term couples the states (b) 
and (d) without creating any defect in the orbital or- 
der. This term therefore enables true coherent motion 
in the insulating ground state with AO order^ Finally, 
the term ^ has yet another effect in that it connects the 
states (b) and (e), as well as the states (e) and (f). In 
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FIG. 1: (Color online) Artist's view of hole propagation in 
the ground state |AO) with AO order presented in (a), with 
vertical (horizontal) bars standing for a (b) occupied orbitals. 

(b) A hole doped into the physical vacuum |AO) at site i re- 
moves four exchange bonds. It can be shifted to new positions 
by consecutively applying different terms of the Hamiltonian: 

(c) a single NN hopping t, (d) a single three-site term r ([8}, 
(e) a single three-site term r and (f) several steps which 
create more flipped orbitals; for more details see text. Broken 
bonds that cost the energy ^ J each are marked by x . 



other words this term connects string states whose num- 
ber of defects differs by two. In some special cases the 
term ([8| may 'split off' clusters of misaligned orbitals. 
Consider for example the state (f) and assume that the 
hole moves upward to site i by virtue of hopping r given 
by Eq. ([8]). Then a 2 x 2 cluster of misaligned orbitals, 
all inverted with respect to the AO order, remains next 
to the hole in the final state. 

In order to discuss the hole motion we restrict the 
Hilbcrt space to a basis of string states^ - — which are 
created by successive application of Tit starting from the 
state (b). Unlike the case of a quantum antiferromagnet, 
the Hamiltonian ((5]) does not produce quantum fluctu- 
ations of the orbital order so that this restriction rep- 
resents an even better approximation than in the spin 
t-J model. Due to the directional nature of Ht and the 
orbital Neel order each hop along the string must be per- 
pendicular to the preceding one so that the maximum 
number of different strings created after n hops is 2". 
The actual number of topologically different strings is 
less than this because one has to exclude self-intersecting 
paths. If we denote the sequence of sites visited by the 
hole by T = {io, ii, ■ . ■ i n } and introduce the 'orbital-flip 
operator' at site i, 



o(i); 



(11) 



where o{i) {o(i)) denotes the orbital at site i which is 
occupied (unoccupied) in the reference |AO) state, the 
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FIG. 2: Numerical results obtained for 20 bands with lowest 
energies, formed by propagating eigenstates |$g (a) en- 
ergy dispersions along the (0,1) direction in the 2D Brillouin 
zone, and (b) the spectral function (|15[) for representative 
values of kb- Parameters: J = OAt, t = O.lt, S = 0.02t. 

corresponding string state can be written as 

\*i,r)=c inAin) J] S,|AO), (12) 

where it is understood that i = io and J- = T' U {i n }- 
Since we want to study coherent hole motion, we con- 
struct Bloch states out of string states (fT2"|) . 

= £ <' k " In . (13) 

3 



Fig. [U^a) shows the dispersion of the lowest bands, 
eh >n (k) as a function of where k = (0, fcf,) is the vec- 
tor from the 2D Brillouin zone. The bands are obtained 
for hole doping at B sublattice and for representative 
parameters^ J = OAt, r = OAt. The energy of the AO 
state with no holes, shown in Fig. [Ha), has been used as 
the reference energy. In agreement with the ID nature 
of the three-site effective hopping ((HJ , the dispersion of 
the lowest energy state and of some excited states shows 
only dependence on kb- 

Figure [21b) shows the single-particle spectral function, 
defined as 

Afru,) = --3^ Zw £Lvi » ^ 

7T ' U) — n (K) + 10 

Z n (k) = |a^ 0jk | 2 , (16) 

where the string Fq corresponds to the state (b) in Fig. [TJ 
In other words, this is just the weight of the bare hole in 
the wave function (fT4"|). As expected, the dispersionless 
bands in Fig. [2fa) have practically no spectral weight 

- as will be seen below, however, these bands give a 
dominant contribution to the optical conductivity. 

The spectral weight of the lowest peak — which would 
form the quasiparticle band at finite doping — shows a 
weak k-dependence which can be understood as follows: 
For momenta near the minimum (maximum) of the dis- 
persion, the hole gains (looses) energy by propagation. 
Since the dominant mechanism of propagation is the hop- 
ping of the bare hole via the three-site hopping term (|SJ) 

— see the transition between states in Fig. [T]from (b) to 
(d) — this gain (loss) in energy will be larger if the weight 
of the bare hole in the wave function is larger. The weight 
of the bare hole, however, also gives the spectral weight 
of the quasiparticle peak. Therefore the weight of the 
peak is larger (smaller) near the minimum (maximum) 
of the dispersion. 



Here we have introduced an additional sublattice index 
S E {A,B}, whereby it is understood that i £ S and 
the sum over j extends all translations of one sublattice. 
The coefficients {ctg 1 ^ k } are variational parameters and 
n stands for a band index. An analogous ansatz for the 
related problem of a hole in a quantum antiferromagnet 
was used before in Refs. I36U381 . In practice, all different 
sets T up to a maximum number n max of defects are 
generated by computer and the Hamiltonian matrix is 
set up. 

We have used n max = 10 and verified that the results 
for the low energy bands are well converged with respect 
to n max . The matrix elements of the Hamiltonian for the 
pairs of states {(b), (c)}, {(b), (d)} and {(b), (e)} are t, r, 
and — r respectively. The sign change with respect to 
the Hamiltonian (|5|) follows from the transformation of 
hopping terms to the hole picture. Due to the directional 
nature of the hopping terms the resulting band structure 
depends on the sublattice index S. 



IV. OPTICAL CONDUCTIVITY 

A. Numerical analysis 

In the next step we discuss the optical conductivity, 
which is defined as a a (uj) = a^aito) + (JB,a(w), with the 
conductivity for sublattice S € {.4, B}: 



O) 



k n=l 

x n sk 5{u- (esn(k) 



(0) 
5,k 



e 5 o(k))}. (17) 



Here a £ {a, b} denotes the direction of the current op- 
erator j a , l^^k) are the approximate single- hole eigen- 
states (QjJ and es„(k) are the corresponding eigenvalues. 
The nsk denote the ground state occupation numbers of 
these states which we assume to be different from zero 
only for the lowest band labeled by n — 0. For a given 
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level of hole doping x they are determined by adjust- 
ing the Fermi energy. This implies that we are assuming 
that for finite hole concentration the lowest band for each 
sublattice is filled according to the Pauli principle. This 
procedure is reasonable for low density of doped holes x. 

We proceed to the discussion of the current operator 
j a (for more clarity the index a is skipped below). For 
the original Hamiltonian (fTJ) this is given by 

J = ** £ ( 6 *> C l+5b,a C ^ + 5k C i+6&,b C i,b) ■ ( 18 ) 
5=±1 

At this point one has to bear in mind that the wave func- 
tions (fT4")l are (approximate) eigenstates of the strong 
coupling Hamiltonian ([5]) rather than the original model 
Eq. (TT]). It is well known that in order to obtain consis- 
tent results for a system described by the original Hamil- 
tonian ((T|) it is necessary to subject the operator in ques- 
tion — here the current operator — to the same canon- 
ical transformation (0| as the Hamiltonian itself. This 
property has been pointed out in the strong coupling ex- 
pansion for the spin Hubbard model^ and we follow here 
this procedure. 

The result of a similar calculation for the present or- 
bital problem is the strong coupling current operator: 

Jt 2g = v(j t +jV+ji?)v + o(J^), (19) 

jt = -«E( a i^ + 6 i c i+b, fl - Ec ')( 20 ) 

i 

i 

- 2 b-E( C U^ +b ,a- H -)' (21) 

i 

isf = =F( a - b ) iT ( C ]±^ a C i,a, C \, 6 C i±a,6 " H -C-) 

i 

T (& + b)ir £ (e^/ i)a ct b c i±&>b - H.c.) . (22) 

i 

It is instructive to trace back the origin of nontrivial 
terms appearing in Eq. (|19l) . The terms which are of 
order t 2 /U stem from processes during which double oc- 
cupancies are created at the intermediate stage due to 
the NN hopping in Eq. (TT5)) and later removed by terms 
in the operator U which are of order t/U or vice versa. It 
turns out that the process during which an electron with 
b orbital flavor moves from one site to its neighbor and 
back, see sites i and j in Fig. [IJa) , does not contribute to 
the transformed current operator j t2 due to the cance- 
lation which originates from the sign dependence of the 
prefactor in Eq. (TT5)) related to the hopping direction. 
Since there is no such dependence in the case of hopping 
term in the Hamiltonian (fT]), the mentioned process gives 
rise to the exchange term in the effective Hamiltonian. 

The above cancelation does not take place in the pro- 
cesses with transfer the hole by two lattice spacings. In 
the first process b orbital moves horizontally, as from site 
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FIG. 3: (Color online) Contribution to the optical conduc- 
tivity <7bs(u) measured along the 8 = a (b) direction (solid 
line) and along the b (a) direction (dashed line), as obtained 
from transitions between states propagating along the a (b) 
direction at hole doping x = 0.1. The total spectrum involv- 
ing transitions from both kind of states is the sum of both 
contributions. Parameters: J — OAt, r = O.lt. Lorentzian 
broadening of width O.Olt has been used. 

m to site I in Fig. [ljb), followed by another hop in the 
same direction onto an empty (hole) site — this gives 
rise to the contribution (|21[) to the effective current. In a 
second process, another b orbital moves horizontally, and 
another hole is created along the diagonal of a plaquette 
(not shown). Such processes give rise to the contribution 
(f2"2")l to the effective current. Now, in the same way as for 
terms in the current operator related with NN hopping 
pop , we may deduce the matrix element of the current 
operator for terms related with further hopping in the 
direction a (or b). For the pairs of states shown in Figs. 
DIb), Did) and in Figs. D^b), Ele), one finds the matrix 
elements — 2ir (0) and ir (—it), respectively. 

The transformed current operator (|19p is next used to 
compute matrix elements between string states. For ex- 
ample, we find that the matrix element for the b compo- 
nent of the current operator jb (|19p between states shown 
in Fig. \V[b) and [lie) is it. The hole shift occurs down- 
wards. The overall sign of the matrix term is positive, 
as the string states (b) Ci. b|AO), and (c) Cj a c\ _ a c i b |AO), 
are defined in the hole language which brings about an 
additional sign change with respect to the prefactor ap- 
pearing in Eq. (1181) . The matrix element of j a for the 
same pair of states vanishes. 

The spectra of optical conductivity (fT7| (obtained by 
applying the Lorentzian broadening width of O.Oli) are 
presented in Fig. [3] The solid line depicts the contribu- 
tion to the conductivity measured along the S — a direc- 
tion from transitions between states propagating along 
the same direction, while the dashed line depicts the con- 
ductivity measured along the 8 — b direction. The true 
response is the mixture of contributions from states prop- 
agating in both directions and for both 5 = a and 6 = b. 
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It is given by the superposition of the spectra plotted 
spectra in Fig. [31 



B. Physical picture of transitions 

Next we give a brief discussion of the physical signifi- 
cance of the optical transitions. In the discussion of hole 
motion we have seen that the dominant hopping term 
T-Lt does not lead to the coherent propagation of a hole 
because it creates a string of defects in the orbital order 
whence the energy increases linearly with the number of 
hops, i.e., with the string length. Coherent propagation 
is enabled only by the conditional NNN hopping term 
(|8|). Let us assume for the moment that this term is 
switched off. Then we can think of localized eigenstates 
of the remaining Hamiltonian 

hM = £ «£W>- (23) 

T 

There are several symmetry operations which transform 
the states \^%,v) into one another: inversion, rotation 
by 7r and reflection by the b (a) axis for i £ A (i £ 
B). This corresponds to the symmetry group C% v and 
the local eigenstates (1231) accordingly realize irreducible 
representations of this group. Next, the Bloch states ([LT| 
may alternatively be written as follows, 

*S)=£ c £ l) l*M>. ( 24 ) 
using the short-hand notation, 

l^) = y|£ <' k ' ! In «l>„ . (25) 

j 

This formulation — which is completely analogous to 
a local combination of atomic orbitals (LCAO) ansatz 
comprising s-like, p-like, d-like basis functions, etcetera 

- immediately clarifies the nature of the optical tran- 
sitions: these are simply dipole-like transitions between 
the approximate local eigenstates \^i, v ) generated by the 
interplay of hopping term H t and the 'string potential'. 
Since the current operator is e.g. odd under rotation 
by 7r, a matrix element {^i^\j a \^i,i>) is different from 
zero only if the two states, l^i,^) and \^fi. v ) have oppo- 
site parity. Since, moreover, the band with the lowest 
energy, n — in Eq. (|14l) . has the totally symmetric 
ground state of the local Hamiltonian as its largest com- 
ponent, the peaks in the optical spectra give essentially 

- with only a small broadening due to weak dispersion 

- by the excitation energies of the states with odd par- 
ity. A very similar interpretation was also given^ i 38 ' 39 for 
the 'mid infrared' spectral weight observed in numerical 
studies of the spin t-J model 3 and applies possibly to 
cuprate superconductors. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the optical conductivity in the 
2D orbital model as a generic model for studying orbital 
polarons in doped orbitally ordered (t2 g or p) systems and 
capturing the essential physics. That effectively spinless 
model is applicable to planar systems like some transition 
metal oxides when (due to the tetragonal crystal field) 
degenerate yz and xz orbitals are active and singly occu- 
pied in a ferromagnetically ordered plane, or to cold atom 
systems in optical lattices whenp orbitals are active. It is 
demonstrated that in the presence of strong electron cor- 
relations, the tendency towards confinement determines 
the properties of the model, i.e., it can be viewed as a sys- 
tem of weakly coupled potential wells. The potential well 
physics originates from exchange energy increase induced 
by the sequences of orbitals flipped by a hole introduced 
by doping, when it moves in the orbitally ordered Mott 
insulator and generates a string potential along its path. 

We have shown that propagating bands are much nar- 
rower in the present orbitally ordered system than in the 
spin t-J model. Results of the analysis based on a variant 
of exact diagonalization, motivated by the string picture, 
suggest the formation of bands in the optical spectrum 
within the Mott-Hubbard gap, by analogy with the mid- 
infrared band observed in doped cuprates^ We predict 
that the optical conductivity would have this form in 
weakly doped Mott insulators with active orbital degrees 
of freedom, while at higher doping orbital stripes would 
form4S Similar to the spin dynamics of stripes in super- 
conducting cuprates^i one expects qualitative changes 
in the optical conductivity for systems with domains of 
AO order separated by orbital stripes, which is an in- 
teresting topic for future studies. Other challenges are 
posed by orbital superfluidity in the p-band of a bipar- 
tite optical square lattice investigated recently^, or by 
spin-orbital systems, where an orbiton may separate from 
a spinon and propagate through a lattice as a distinct 
quasiparticle.— 
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Appendix: Origin of optical transitions 

The aim of this section is to analyze the mechanism 
underlying the doping induced formation of states lying 
within the Mott-Hubbard gap of a correlated insulator 
which is orbitally ordered. Furthermore, it will be dis- 
cussed, how the structure of those states influences the 
optical response. We believe that despite some simplifica- 
tions, conclusions which will be drawn from this analysis 
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FIG. 4: (Color online) The mechanism of the Trugman 
process^ allowing for hole deconfinement in a doped anti- 
ferromagnet described by the Ising model, (a) A hole doped 
to the Mott insulator with the antiferromagnetic order (spins 
are represented by arrows) generates six bonds with parallel 
spins by interchanging its position with spins on a plaquette 
after three hops. By a clockwise/anticlockwise hopping by 
three more steps one arrives at states (c)/(b), where all spin 
excitations are removed and the antiferromagnetic state is re- 
paired — hence such processes couple the configurations (b) 
and (c) . Broken bonds that cost the classical energy | J each 
are marked by x. 



are applicable to more complex situations encountered in 
real systems revealing orbital order. 

Due to the simplicity of the analyzed model the rele- 
vant part of the Hilbert space consists of states which can 
be obtained by subsequently shifting the hole created in 
the AO state, the situation shown in Fig. [ljb), to nearby 
sites. Thus we could use that original position of the hole 

— R(R, F) [here R(R, F) stands for the function of the 
present hole position R and the set of sites on which or- 
bitals have been flipped (J 7 )] to label a given string state 

^E'r.j 1 '}- In order to determine R(R, F) in a unique way, 
we additionally demand that the total length of hole path 
necessary to create the state I^r^) from a state repre- 
senting a hole created at the site R(R, F) in the AO state 
is minimal. For example, in the case of the state shown 
in Fig. [ljf) the vector R refers to the site n while the 
vector R(R, F) to the site i. In principle there could be 
more than one "original site" obeying those conditions 
for a given pair {R, F}. In that case we would arbitrar- 
ily choose a single R(R, F) and the new method of state 
labeling by the new pair {R(R, F), F} would also work. 

The possible ambiguity in determining the original po- 
sition of the hole for a given string state would give rise 
to a new channel of coherent hole propagation which is 
known in the case of doped antiferromagnets. Fig. [4] de- 
picts the so-called Trugman process^ which brings about 
hole deconfinement on a square lattice even in systems 
with Ising-type anisotropic exchange interaction as in the 
present orbital t- J z model. The state shown in Fig. QJa) 
can be obtained both by shifting the hole clockwise or an- 
ticlockwise around the elementary plaquette by three lat- 
tice spacings (initially the hole replaced a 4,-spin as shown 
in Figs. SJb) and|4jc), respectively). In other words, by 
performing one and a half of the circular movement of the 
hole on a given plaquette, the end effect is that the hole 
has moved along the plaquette diagonal without bringing 



about spin flips in the Neel state, which gives rise to the 
weak coherent hole propagation. 

Now, with the help of the computer algebra, we will 
demonstrate that in the present two-orbital problem the 
ambiguity in determining the original position of the hole 
for a given string state does not exist in the case of a 
string state generated by a single hole, with a small num- 
ber of twisted orbitals in the |AO) state. Thus, it seems 
that the only channel allowing for hole propagation is due 
to the Hamiltonian induced coupling between states rep- 
resenting a hole created in the |AO) state, as the states 
depicted in Figs, [ljb) andQJd). This kind of hopping is 
mediated by the small r -c t term flSJ. 

We proceed now to provide the justification of that 
statement when restricted to the low energy sector of the 
Hilbert space, i.e., consisting of states with limited num- 
ber of flipped orbitals in the AO) state. The Hamiltonian 
matrix represented in terms of states l^s^k) consists of 
two decoupled blocks. The first (second) block is formed 
by states coupled with the state representing a hole cre- 
ated in the perfect |AO) state at a site belonging to the 
sublattice A (B). Despite that both blocks show explicit 
dependence on both {k a , kb}, the energies of eigenstates 
for each block disperse along a single direction, either 
(10) or (01) in the 2D Brillouin zone. 

By representing the Hamiltonian matrix in terms of 

|*s^,k> - M E^^l**^)' (A.l) 

which is equivalent with performing a kind of gauge 
transformation, we get rid of the superfluous momentum 
dependence of the Hamiltonian matrix. Furthermore, if 
we neglect the term (jHJ, the Hamiltonian matrix lacks 
any dependence on momentum, which shows that the 
hole deconfinement occurs here solely due to that term 
and that the hole becomes confined when r is set to zero. 
Those findings are restricted to the basis considered by 
us which is limited by the path length. Furthermore, we 
do not probe the part of the Hilbert space which is not 
coupled by a Hamiltonian power to original states, repre- 
senting a hole created in the |AO) state. Thus, a general 
mathematical proof of the hypothesis regarding the na- 
ture of the propagation and valid for the whole Hilbert 
space is still needed. 

To identify the origin of optical transitions contribut- 
ing to the spectrum depicted in Fig. [3l we analyze now 
the properties of the Hamiltonian matrix represented in 
terms of states (jA.lj) [i.e., in the case when the original 
hole position is used to label string states], when r is 
set to 0, which means that only the "fast" part of the 
Hamiltonian is considered. It turns out that this part 
determines the overall structure of the energy-band hier- 
archy. Due to the lack of energy dispersion the bands are 
now completely flat, but their positions correspond very 
well to the sequence of bands obtained for finite r and 
shown in Fig. [2J The picture which emerges from that 
correspondence is that the physics of hole hopping in the 
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FIG. 5: (Color online) Artist's view of some disconnected 
string states which are not coupled by NN hopping t with 
the original state of a hole doped in the AO state shown in 
Fig. []Jb). After moving by one hopping process t from site 
i to site j, the hole hops downwards by the effective hopping 
t and leaves behind an island of excited bonds (a). Further 
hopping of the hole by either (b) one t step, or (c) two t steps, 
generates more defects in the |AO) state. Broken bonds (x) 
as in Fig. 4. 

orbitally ordered background at large energy scale is de- 
termined by fast moves to NN sites, accompanied by the 
creation of defects in the orbital arrangement. Due to 
the increasing length of defect sequences (strings) a hole 
behaves like a particle in a potential well. Consequently, 
band hierarchy corresponds to the sequence of eigenener- 
gies for the corresponding problem of the particle in the 
well. On the other hand, the energy dispersion is deter- 
mined by "slow" hopping which brings about the modifi- 
cation at the energy scale ~ r <C t and can be viewed as 
a perturbation introduced on top of the robust structure 
of eigenstates arising for the potential-well problem. For 
finite r, the Hamiltonian matrix block formed by states 
coupled with the state representing a hole removed from 
the A (B) sublattice in the |AO) state does not show 
any dependence on k a {kb), which explicitly demonstrates 
that the hole propagation is ID. 

The point group of the underlying orbital background 
is C^v The Brillouin zone gets folded due to the stag- 
gered form of the AO order and all parameters like the en- 
ergy are periodic with the periodicity (ir, w) and (tt, —it). 
The bottom of the lowest energy band lies on the lines 
(7r/2, kb) and (fc ,7r/2) (and on lines equivalent by sym- 
metry) for two orthogonal directions of the ID hole prop- 
agation, respectively. Within the single particle approx- 



imation which we apply here, the Mott insulator doped 
with holes starts to fill orbital polaron bands near their 
bottoms. Even for nonzero r the Hamiltonian blocks rep- 
resented in the basis of states (|A.1[) are fully symmetric 
with respect to the point group C% v for wave vectors lying 
at the band bottom. For example, the block in the sector 
consisting of states lacking dispersion in the b direction 
is symmetric at the wave vector (7r/2, kb) with respect to 
the reflection in the a axis, as it lacks the dependence 
on kb- The inversion in the b axis transforms (7r/2, kb) 
into (— 7r/2, kb), which is the same as (7r/2, kb + n) due to 
Brillouin zone folding, and equivalent to (n/2, —kb) due 
to the lack of dependence on kb- Thus we can use Ci v 
to classify the symmetry properties of states at the band 
bottom4^ 

The ground state for both nonzero r (system with 
propagating holes) and for t set to zero (system with 
confined holes) is even with respect to both reflections: 
in the a axis (s a ), and in the b axis (s&). In the latter case 
the states like the ones depicted in Figs. [5][a)-[5][c) which 
are not coupled by a product of the t-hopping terms Q 
with the original state, with a hole created in the |AO) 
state at site i and shown in Fig. [Hb), do not contribute 
to the ground state, while in the former case their weight 
is small. The state depicted in Fig. EJa) is coupled with 
the state shown in Fig. [TJc) by the r-hopping term ((SJ) 
which, on the other hand, is coupled with the original 
state presented in Fig. QJb) by the t- hopping term |6]). 
The states shown in Figs. Etb) and Etc) are both coupled 
by hopping terms with the state presented in Fig. EJa). 

The first excited state in the sector corresponding to 
the propagation in the a direction is odd with respect to 
the reflection s a and even with respect to Sb- Since the 
current operator (p~9|) is an axial vector, its b component 
(which is odd with respect to s a ) couples the first excited 
state with the fully symmetric ground state, which gives 
rise to the contribution to the optical weight in the form 
of the first peak from the left shown in Fig. [3] Upon 
doping, the optical transitions occur at momenta located 
in the vicinity of the band minimum, which explains why 
the peak gets a finite width. In the system with no prop- 
agating holes, the states such as those depicted in Fig. 
El which are not coupled by NN hopping to the original 
state, do not contribute to the first excited state, while 
their weight is small when the holes propagate. 
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